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Abstract — The Station Explorer for X-ray Timing and Navi- 
gation Technology (SEXTANT) is a NASA funded technology- 
demonstration. SEXTANT will, for the first time, demonstrate 
real-time, on-board X-ray Pulsar-based Navigation (XNAV), a 
significant milestone in the quest to establish a GPS-like nav- 
igation capability available throughout our Solar System and 
beyond. 

This paper describes the basic design of the SEXTANT system 
with a focus on core models and algorithms, and the design 
and continued development of the GSFC X-ray Navigation 
Laboratory Testbed (GXLT) with its dynamic pulsar emulation 
capability. We also present early results from GXLT modeling 
of the combined NICER X-ray timing instrument hardware and 
SEXTANT flight software algorithms. 
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1. INTRODUCTION 

SEXTANT is a NASA Space Technology Mission Direc- 
torate Game Changing Development Program Office funded 
technology-demonstration enhancement to the Neutron-star 
Interior Composition Explorer (NICER) mission. NICER 
is a NASA Science Mission Directorate funded X-ray As- 
trophysics Explorer Mission of Opportunity that will be 
hosted on the International Space Station (ISS). SEXTANT 
will, for the first time, demonstrate real-time, on-board X- 
ray Pulsar-based Navigation (XNAV), a significant milestone 
in the quest to establish a GPS -like navigation capability 
available throughout our Solar System and beyond. The 
SEXTANT demonstration will exploit the large collecting 
area (> 1800 cm 2 ), low background (< 0.2 counts/s), and 
precise timing (< 100 ns 1-cr) provided by NICER’s X-ray 
timing instrument. 

SEXTANT will demonstrate the use of Millisecond Pulsars 
(MSPs), rapidly spinning neutron stars, as deep-space navi- 
gation beacons which could someday guide humankind out 
of Earth orbit, to destinations throughout the Solar System 
and beyond. SEXTANT’S primary objective is to demonstrate 
real-time orbit determination with uncertainty better than 
10 km with 2 weeks of measurements in the highly dynamic 
Low-Earth Orbit (LEO) of the ISS. SEXTANT also includes 
the development of a ground testbed, called the GXLT, that 
enables real-time and faster than real-time simulation of nav- 
igation scenarios. With the addition of a unique Modulated 
X-Ray Source (MXS), flight-like X-ray detector, and time- 
tagging electronics, this testbed provides a test-as-you-fly 
hardware-in-the-loop simulation capability. 

The SEXTANT demonstration is a technology enhancement 
to the NICER mission, which is an X-ray Astrophysics 
Mission of Opportunity to ISS that is currently in Phase 
C [1] and plans to launch in October 2016. NICER will 
undertake a fundamental investigation of extremes in grav- 
ity, material density, and electromagnetic fields of rapidly 
spinning neutron stars via time-resolved X-ray spectroscopy 
[2, 3]. NICER achieves this objective by deploying an X-ray 
telescope instrument as an attached payload on a zenith- side 
ExPRESS Logistics Carrier (ELC) aboard the ISS, Figure 1. 
NICER offers over an order-of-magnitude improvement in 
time-coherent sensitivity and timing resolution beyond the 
capabilities of any X-ray observatory flown to date. As a 
flight software augmentation to NICER, SEXTANT will use 
a subset of the data collected for the NICER science program, 
on-board and in real-time, to accomplish its objective. 

NICER and SEXTANT are funded through a cost-sharing 
opportunity between the NASA Science Mission Directorate 
and NASA Space Technology Mission Directorate Game 
Changing Development Program Office. 

In the following, an introduction to X-ray pulsar navigation 
is given, followed by a brief description of the NICER X- 
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Figure 1: Neutron- star Interior Composition Explorer 
(NICER) payload reference representation. 


ray timing instrument and the SEXTANT architecture. Then, 
details of the SEXTANT core models and algorithms, and 
end-to-end simulation are discussed. This is followed by a 
description of the GXLT hardware components and a brief 
presentation of prelimiary simulation results. Finally, the 
paper concludes with a summary of ongoing and future 
work. 


2. X-Ray Pulsar Navigation 

X-ray observations of celestial sources can provide useful 
navigation information to spacecraft in a range of applications 
from LEO to interplanetary, and even interstellar, space. 
One source of such information are X-ray emitting pulsars, 
which are neutron stars whose X-ray emission is modulated 
at the rotational period of the star. In this paper, we define 
the X-ray Pulsar Navigation (XNAV) concept as spacecraft 
navigation using X-ray observations of such pulsars. A subset 
of pulsars, the millisecond pulsars, are highly stable clocks, 
with long term stability comparable to laboratory atomic 
clocks. For these pulsars, a simple physical model with a 
small number of parameters can predict the arrival time of 
pulses to microsecond accuracy over months or years. A 
measurement of the difference between the arrival time of a 
pulse at a spacecraft and the predicted arrival time according 
to an onboard navigation solution can provide an error signal 
that can be used to measure the location of the spacecraft in a 
manner similar to GPS [4-6] (see Figure 2). 

The use of radio pulsars as navigation beacons was first 
considered shortly after their discovery [7]. The idea was 
later extended to X-rays using the earliest established X-ray 
pulsars [8], but the achievable accuracy was severely limited 
by the noise characteristics of the X-ray pulsars known at 
the time. The first X-ray instrument with a specific goal of 
exploring X-ray navigation techniques was the Unconven- 
tional Stellar Aspect (USA) Experiment, flown in 1999 on 
the DoD ARGOS satellite, under the Space Test Program 
[9]. This experiment explored a broader concept of X-ray 
navigation, not limited to pulsars and also not limited to 
position determination but considering also time transfer and 
aspect determination. 


During the 1990s X-rays began to be detected from millisec- 
ond pulsars previously known only as radio pulsars [10]. This 
development greatly improved the expected accuracy of an 
XNAV system and spurred detailed studies [4], resulting in a 
patent on the idea (US Patent 7,197,381). A DARPA program 
emphasizing millisecond pulsar methodologies represented 
the next stage of DoD development, and it was during this 
program that the first laboratory facility to simulate X-ray 
pulsars was developed at GSFC. By now the concept of X- 
ray navigation is being pursued in several other countries as 
well as the U.S [11, 12]. 

XNAV has the potential to become an enabling technology 
for very deep space exploration missions and an important 
augmentation to NASA’s Deep Space Network (DSN), the 
current standard for interplanetary navigation and communi- 
cation [13]. 


3. The NICER X-ray Timing Instrument 

The key measurable for an XNAV instrument are pulse arrival 
times determined from a set of detected X-ray photons. The 
instrument must be designed to be able to produce a high 
signal-to-noise pulse profile in an integration time that is 
short compared to the timescale for a propagated orbit to 
deviate from the true trajectory by more than the required 
navigation accuracy. In turn, this requires a detector with the 
following characteristics: large effective area in the region of 
the X-ray spectrum where MSP pulsations can be observed 
(roughly 0.2-8 keV), high precision time tagging of each 
X-ray photon, and low background rates. To be useful 
in future spacecraft navigation applications, this should be 
achieved with the lowest possible mass, volume, and power 
requirements. 

The NICER X-ray Timing Instrument (XTI) (Figure 3) is 
extremely well suited to this task. It is a modular array 
of 56 identical telescopes, making it easily scalable to a 
range of potential applications. Each telescope consists of 
a lightweight grazing-incidence optic made up of concentric 
foil mirrors. The mirrors concentrate X-rays onto a small (1 
mm radius unobstructed circular aperture) detector area using 
a single bounce, in contrast to typical (e.g. Wolter I; [14]) 
imaging X-ray optics that require two bounces and thus incur 
a significant efficiency penalty to achieve quality imaging. 
The XTI provides over 1800 cm 2 of area in a package with 
frontal area of 6400 cm 2 , an areal efficiency of 28%. 

The X-ray detectors are commercial (Amptek) silicon drift 
detectors that have very high quantum efficiency over the 
photon energy range of interest. The detectors are read out 
by dual-channel electronics chains that provide both high 
time resolution (100 ns) and excellent spectral resolution (120 
eV) with very low dead time. Particle backgrounds are low 
because of the small detector volume made possible by the 
concentrating optics, while particles that do interact in the 
detector can be rejected at high efficiency by their energy 
deposition and by filtering events that occurred outside of 
the illuminated area of the detector by a comparing the pulse 
heights determined by the slow and fast electronics channels. 
In addition to reducing particle backgrounds, the focusing 
of the mirrors reduces the backgrounds from the cosmic 
diffuse X-ray background and neighboring sources that are 
outside the ~ 6 arcmin (FWHM) field of view. The expected 
background rate in the critical 0.4-2 keV band is < 0.2 counts 
per second. 
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Figure 2: Geometry of an XNAY observation. 



Figure 3: Block diagram of NICER XTI showing the main 
components: concentrators, focal plane modules, detectors 
and measurement and power unit. 


4. SEXTANT System Architecture 
Overview 

The main components of the SEXTANT system and the 
data flow relationship among these components are shown in 
Figure 4. A detailed overview of the SEXTANT architecture 
can be found in [15]. The components that are external to the 
NICER XTI are briefly summarized below. 

SEXTANT flight software 

The X-ray Pulsar Navigation Flight Software (XFSW) is a 
single application hosted by the NICER Instrument Flight 
Software (IFSW). It pre-processes and buffers photon events 
from the XTI until a sufficient number of events from a 
single pulsar is collected. The photon events are then batch 
processed using a Maximum Likelihood (ML) estimation 
algorithm to produce a single pulse phase and Doppler 
measurement. The Goddard Enhanced Onboard Navigation 
System (GEONS) navigation filter [16], seeded by a single 
(intentionally degraded) initial state estimate from the NICER 
GPS receiver, blends models of spacecraft dynamics with 
these pulse phase and Doppler measurements to maintain 
an estimate of the relevant spacecraft state, i.e., position, 
velocity, drag coefficient, etc. 

SEXTANT ground system 

The SEXTANT ground system maintains the pulsar catalog 
that provides current timing models, or ephemerides, and 
pulse profile templates as needed for navigation processing. 
The ground system is also responsible for performance moni- 
toring and telemetry collection for post-processing purposes. 

The pulsar catalog consists of a list of SEXTANT pulsars, 
with their associated timing models and X-ray lightcurve tem- 
plates. The current baseline SEXTANT pulsar list is provided 
in Table 1, along with some basic parameters. Templates for 


six of the catalog pulsars, as constructed from archival X- 
ray data retrieved from NASA’s High Energy Astrophysics 
Science Archive Research Center (HEASARC) [17], are 
shown in Figure 5. The ground system will measure the phase 
relationship between the radio and X-ray templates for each 
pulsar, and track variations in the pulsar dispersion measure 
that are needed to maintain this alignment. 

The timing models will be regularly updated using Time-of- 
Arrival (TOA) data from external observatories, and com- 
bined with XTI data once NICER is operational. The ground 
system relies on the Tempo2 pulsar timing software [18, 19] 
for the generation of timing models by fitting parameterized 
models to measured radio and X-ray pulse TOAs, and for 
generating the polynomial fits to the full models that are used 
for onboard processing. 

SEXTANT ground testbed 

The GXLT is a unique software and Hardware-in-the-Loop 
(HWIL) environment for development and testing of various 
XNAV algorithms and detector technologies. This allows the 
main SEXTANT core models and algorithms, as described in 
Section 5, to be exercised in a flight-like environment without 
a full-scale NICER XTI. The ground testbed is described in 
detail in Sections 7 and 8. 


5. SEXTANT CORE MODELS AND 
Algorithms 


Photon arrival model 

The fundamental data provided to the SEXTANT XNAV 
algorithms are the observed photon event timestamps (in GPS 
time, accurate to 100ns RMS) with associated pulse heights, 
proportional to photon energy. In this section, we present the 
formal statistical model used to model this process. 

The observed photon event timestamps {T^}^ =1 received 
on-board NICER from the XTI are modeled as the arrival 
times of a Non-Homogeneous Poisson process (NHPP) with 
time varying mean cumulative count function A (£) [20, 21]. 
For such process, the probability of k events arriving in the 
interval (£ a , tf) is given by 


P{k\ {t a ,t b )) 


e -(A(t 6 )-A(0)(A(t 6 ) - A(ta)) k 
~k\ 


( 1 ) 


We assume the existence of a “rate function” A such that 
A (£) = f Q A (s)ds. We further assume that the rate function 


3 



tit t lt.lt 


GPS time stamped 
Photon + Background 
events, GPS telemetry 


Ground System 
Software 


Ground System 



GPS time stamped 
Photon + Background 
events 


Flight Software 


Event 

Filtering 


Pulse 

Phase 

Estimation 


GEONS 

navigation 

processing 


Pulsar almanac entries (timing models, 
templates) and alert messages 


1 ±_ 

Pulsar almanac j 

TLM 

* 

► 



Test stimulus 


Lab Testbed 

Goddard XNAV 
Lab Testbed 
(GXLT) 


Test stimulus 


Pulsar data from observatories 
(Arecibo, FERMI, etc.) 


Proving ground for algorithm 
development, validation, test 


Figure 4: SEXTANT system architecture. 

Table 1: SEXTANT pulsar catalog with basic parameters 


Name 

Period ( P ) 
(ms) 

D a 

(kpc) 

Source Pulsed Rate 
(a, cnts/s) 

Total Bkg Rate 

(/3, cnts/s) 

i P 

(-) 

<j c 

(Ats) 

Crab Pulsar 

33.51 

2.0 

660.000 

13 860.20 

56 963.3 

6.6(10) 6 

B1937+21 

1.56 

5.0 

0.029 

0.24 

27.0 

14.1 

B1821-24 

3.05 

4.9 

0.093 

0.22 

231.2 

9.5 

J02 18+4232 

2.32 

5.8 

0.082 

0.20 

5.7 

45.8 

J0030+045 1 

4.87 

0.3 

0.193 

0.20 

5.4 

98.8 

J1012+5307 

5.26 

0.7 

0.046 

0.20 

— 


J0437-4715 

5.76 

0.2 

0.283 

0.62 

2.9 

159.1 

J2124-3358 

4.93 

0.3 

0.074 

0.20 

— 

— 

J22 14+3000 

3.12 

1.5 

0.029 

0.26 

— 

— 

J0751+1807 

3.48 

0.4 

0.025 

0.22 

— 

— 

J 1024-07 19 

5.16 

0.5 

0.015 

0.20 

- 

- 


a D is the distance to the pulsar in kpc (1 pc = 3.26 light-years) 

b The measurement accuracy for this observation length is 1 /as, but the intrinsic rotational instability 
of the Crab Pulsar imposes a limit on the model prediction accuracy of about 10 /as over a few days. 
c The phase error is given by a — V— with T — 1800 s (obtained from the CRLB, Eq. (9)). 


has the form 

X(t) = \(<f>(t)) = p + ah(<j>(t)) (2) 

with a > 0, P > 0 constant, and h a nonnegative, continu- 
ous and piecewise differentiable real function, periodic with 
period 1 (h((j) + m) = h{^>) for all integers m) and with unit 

area f* h{(j))d(j) = 1 (see Fig. 5 for six pulsar profiles from 
the SEXTANT catalog). Finally, 0 is assumed smooth with 
bounded positive derivative, bounded away from zero, so that 
</> _1 exists and is smooth. 

The parameters ck, /3 represent the average signal and total 
background count rates in units of counts per second, respec- 
tively. The function h describes the characteristic shape of the 
pulsar intensity profile or lightcurve , see Figure 5. Finally, </> 
represents the evolution of the phase of pulsation with respect 


to the time parameter t as seen at the detector. 

Pulse arrival model 

The phase at the detector 0 is modeled as 

0(f) =00 (3) 

where (j ) o is the phase evolution at a hypothetical reference 
observatory and r(t) is the light propagation time of the pulse 
wavefront moving from the detector to the reference observa- 
tory. SEXTANT’S flight software and end-to-end simulation 
currently uses either a Geocenter or Solar System B ary center 
(SSB) reference observatory. In the former case, cr(t) ~ 
n • x(t), where x(t) gives the detector coordinates in a frame 
centered at the reference observatory and ft is the direction 
to the given pulsar, to a sufficiently accurate approximation. 
In the latter case, parallax and Solar Shapiro delay terms are 
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Figure 5: Current set of template profiles for XNAV pulsars. Two cycles of each pulsar are shown. 


additionally needed to achieve similar accuracy. 

SEXTANT relies on the pulsar timing software Tempo 2 
[18] to provide a model for phase evolution at a reference 
observatory. Tempo 2 models are physically motivated least- 
squares fits to radio observatory data of the form 

<t>o (t) =p(t- A(t)), (4) 

with p a quadratic or cubic polynomial shifted by timing 
correction A (t) which includes terms for Roemer, Einstein, 
and Shapiro, and binary delays. 2 Tempo2 also provides 
convenient piecewise polynomial approximations to the full 
timing model 0o, which are used by the SEXTANT flight 
software. 

Count rate models 

The pulsed count rate a observed from each pulsar with the 
NICER XTI is determined by the spectrum of the source 
and the response matrix of the instrument. These values 
will be determined with high accuracy during pre-launch and 
on-orbit calibration for each target in the SEXTANT pulsar 
catalog (Table 1). 

The background count rate is composed of any unpulsed flux 
from the pulsar, other point sources in the field of view, 
diffuse X-ray flux in the aperture, and unrejected radiation 
background. The latter component is expected to vary with 
position and time, and is the sum of the particle (proton, 
electron, positron) interactions and those originating with 
high-energy (> 30keV) X-rays and 7-rays that penetrate 
the instrument enclosure and produce secondary radiation 
within the XTI. The current SEXTANT simulations assume a 
constant background rate over each observation interval, but 
work is underway to remove this simplification and update (3 
to account for on-orbit particle background variability. 


2 The latter term accounts for the orbital motion and relativistic effects for 
pulsars in binary systems, which constitutes the majority of pulsars in the 
SEXTANT catalog. This effect prevents (3) from being well approximated 
by a linear function in time over typical observation timescales. 


Pulse phase and Doppler estimation 

The baseline SEXTANT algorithms work by observing sev- 
eral MSPs from the SEXTANT catalog cyclically and when 
sufficient time is accumulated on a particular MSP, the col- 
lected photon events are batch processed to extract pulse 
phase and Doppler estimates, which are finally passed to a 
navigation filter. 

We would like to use standard parameter estimation tech- 
niques to generate the phase and Doppler estimates. In 
general, due to the spacecraft motion and the form of 0 0 
(especially with the binary terms), the phase evolution has 
a complicated form that is difficult to directly parameterize 
with a low-order model. For example, a periodic NHPP 
model, used in some prior work, would be inadequate in most 
cases. To address this, as described in the following, our 
approach assumes the existence of an approximation to the 
phase evolution that captures most of the dynamics and we 
then parameterize a correction to this approximation. 

Specifically, we assume existence of a smooth approximation 
0 such that 0(f) = 0(£) +e(£), so that {70} is NHPP with rate 

A (0(£) + e(£)). Furthermore we assume over the observation 
interval [t a ,tb], e(t) can be fit by a low-order parametric 
model. 

The phase prediction 0 

Using a predicted position x(t), provided by the GEONS 
navigation filter, we expand the phase function as 


= 0 O [t- 


x(t) • ft 


( x(t) • n 5x(t ) • ft 

= <PQ\t 1 


- 0(f) + 00 [t - 


x(t) • n \ Sx(t) • ft 


= 0(£) + e(t), 


(5) 


5 


where x is the spacecraft position vector, n is the unit vector 
direction to the pulsar, c is the speed of light, 8x(t) = x(t) — 
x(t), and e(£) := p o (t — x ^' n ^j 6x d)-™ • see observation 

geometry shown in Figure 2. We assume e(£) is well fit by 
a linear model, i.e., e(t) ~ q + /(£ — t a ) for constant q and 

/• 3 

In summary, the SEXTANT phase and Doppler estimation 
routine proceeds as follows. 

1. Observe arrival times {tk}^ =1 during a fixed interval 

\pa i tb \ • 

2. Determine estimates (g, /) of the parameters (g, /) in 
the model e(£) = q + /(£ — £ a ). 

3. Form phase and Doppler estimates p(t) = 0(£) + q + 
/(£ — £ a ) and 0 (f) = /, respectively. 

Step 2 is clarified below. 

Pulse phase estimation routines 


parameters 6 = (g, /) (e.g., [25]). It is shown in [ 6 , 24] to be 
given by 


CRLB(<9) = — 
I p 


2 -3 

r p r r2 

-3 _6_ 

^2 ^3 


with T = tb — t a , and where constant I v is 


( 9 ) 



wm? 

P + ah(p) 


dp. 


( 10 ) 


For SEXTANT, due to a constraint imposed by the navigation 
filter, we always refer the measurement to the time at end of 
the observation interval 4 and get 


4>(t b ) = <f>(t b ) +q + f{t b -t a ). ( 11 ) 

Subtracting this equation from <j>{t b ) and using the assumed 
form of the prediction error we get that the phase and fre- 
quency estimate errors areg — g + (/ — f)T and (/ — /), 
respectively. Thus using (9), we get y4^ and respec- 
tively, for the phase and frequency error variances. 


As described previously, the photon time-of-arrival process 
{Tk }% =1 observed during the interval [£<*,£0 is assumed to 
follow a NHPP with rate 

Mt-,q,f):=X(4>(t) + q + f(t-t a )), ( 6 ) 

and our goal is to estimate the parameters (g, /). Toward this 
end, it can be shown that the arrival phases admit the joint 
probability density function 


Beyond Poisson statistics and the CRLB, there are additional 
sources of error in the derived state vector estimates that 
should be considered. Addition of the following terms are 
being developed as enhancements to the current baseline 
SEXTANT models and simulations. The first three amount 
to errors in the pulsar almanac parameters po, a , 0 , and h. 
We strive to control these error sources to be smaller than the 
Poisson noise contribution, which partially justifies the use of 
twice the CRLB for the filter measurement noise. 


p(t 1 ,t 2 ,...t N ;q,f)~e A(s)ds II^ = 1 A(f fc ) (7) 

with the constraint t a < t\ < £2 . . . < £at < £ 5 , where N is 
also a random variable here (see e.g., [ 6 , 23, 24]). 

An approximate MLE — The exponent in (7) is the mean total 
count rate, which approximately equals (a + /?)(£& — £ a ), 
and is thus considered independent of the parameters (g, /), 
as long as the observation interval includes many cycles 
of pulsation. This is always the case for SEXTANT MSP 
observations, therefore, this term may be dropped leading to 
the approximate MLE 

N 

(q,f) = argmaxV'logA (<p(T k ) + q + f{T k - t a )). ( 8 ) 
qJ 


Timing noise and glitches — Pulsar timing noise and glitches 
(abrupt shifts in pulse frequency) are both unpredictable pro- 
cesses that limit the accuracy of extrapolated timing models. 
These effects are most pronounced for young pulsars (such 
as the Crab pulsar), but even the very stable millisecond 
pulsars usually have measurable timing noise and occasion- 
ally even suffer small glitches (as seen in the XNAV pulsar 
B 1821— 24). For glitches, the only mitigation is to monitor 
the timing of each of the pulsars and stop using it for naviga- 
tion when a glitch is detected, only returning it to operation 
when a sufficiently stable post-glitch solution is determined. 
For timing noise, the effects can be minimized by updating 
the timing models frequently based on recent observations (to 
limit the time the pulsar has to wander away from the current 
solution), and to use optimal techniques for extrapolating 
pulsar phase predictions in the presence of timing noise, as 
described in [26] . 


Measurement uncertainty 

Here we discuss the components of the error in the phase 
and Doppler estimates. An overall estimate of the level 
of measurement noise is a required input parameter for the 
navigation filter. For SEXTANT we specify a measurement 
noise (lcr) level corresponding to two times the square-root 
of Cramer-Rao Lower Bound (CRLB), described next. 

Poisson noise — Given sufficient observation time, the ML 
described above will achieve the CRLB for model (7) with 


3 As suggested in [22], an approximation better matched to the orbital 
dynamics would be Sx(t) = <&(t,t a )Sx(t a ), where <!>(£, t a ) is the state 
transition matrix provided by the navigation filter. In this case, we would 
estimate a 6-dimensional parameter vector Sx(t a ). We leave this for future 
investigation because the linear model has been shown to work well in 
simulation. 


Uncertainty in radio to X-ray offset — Initially, the timing 
models for XNAV pulsars will be developed from ground- 
based radio observations, but the phase measurements are 
made in the X-ray band. Therefore, the absolute phase offset 
between the X-ray and radio pulse profiles must be known 
precisely. This offset can be measured, through contempo- 
raneous measurements with X-ray and radio telescopes, to 
an accuracy limited by signal-to-noise of each measurement 
and the cross calibration of the reference clocks used by each 
observatory. Additionally, the radio pulse arrival times must 
be corrected for the effects of dispersion. This can be done 
with careful measurements at multiple radio frequencies. 
However, any variations from the measured value introduce 
an error in the phase predictions at infinite frequency (e.g. 


4 Although, referring it to the middle of the interval would improve the 
estimate. 
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X-ray). These variations can come from several sources in- 
cluding the interstellar medium (as the pulsar moves through 
the Galaxy at high velocity), the interplanetary medium, the 
solar wind, and the Earth’s ionosphere. 

Astrometric errors — The measured astrometric parameters 
for the pulsars (position, proper motion, and parallax) are 
used to compute the time offset between the spacecraft and 
the reference observatory. These quantities are measured 
via timing (or interferometric) observations at 1 AU from 
the Sun (i.e. on Earth). As the spacecraft moves to much 
larger distances from the Sun, the accuracy required on these 
parameters increases but this is not a concern to SEXTANT. 

Non-timing model mismatch — Errors in the pulse template 
mismatch and source and background flux rates (if they are 
not estimated, which they will be in the final SEXTANT 
algorithms) will introduce additional error in the estimation 
process. 

Non-ideal behavior of the reference clock — The events pro- 
vided by the NICER XTI are timestamped to GPS time to 
100 ns accuracy. The SEXTANT baseline implementation 
uses the full accuracy of these timestamps, effectively pro- 
viding continuous calibration of its onboard clock. An XNAV 
sensor operating outside of the range of GPS will timestamp 
event times with respect to an on-board clock which will 
drift over time. This error can be compensated either by 
periodically recalibrating via contact with the ground system 
or by enhancing the measurement model to allow estimation 
of clock states, or a combination of both. While our baseline 
approach is idealized, it is sufficiently representative for the 
purposes of the demonstration. Moreover, non-ideal refer- 
ence clock behavior will be investigated on the ground by 
intentionally degrading stored science telemetry. 

Event simulation 

The ability to efficiently simulate NHPP photon events is 
a core function of the GXLT necessary for developing and 
testing the SEXTANT flight software algorithms. 

An important property of Poisson processes, relevant to sim- 
ulation, is that given a standard (with constant rate 1) Poisson 
process {Sk}, the process 70 = A ~ 1 (S k ) is an NHPP with 
mean cumulative count function A (£). 


Thus, an NHPP with rate A (t) = A (0(f)) can be generated 
by Tk = A -1 (S^) where {Sk} is a standard Poisson process 
(with constant rate 1) and 

A (t) = f \{4>{z))dz. (12) 

Jo 


However, a direct inversion of A, through table lookup, for 
example, would be computationally expensive and does not 
take advantage of the periodicity of A. Considering a change 
of variable u = f(z) in (12), we have 


A (t) = 


rHt) 

Jcb(O) 


\{u) 


m H4> HA)) 


-du. 


0(^ra) 


MHO) -■ Ai (HO) 


(13) 


where A (t) = ffl X(u)du and we assume without loss of 
generality that 0(0) = 0. We also assume 0 (£) ~ 0 (£ m ) 


holds approximately over the simulation interval where t m is 
the midpoint of the observation interval, and we define Ai := 

A /Ht m ). 


Now we can efficiently simulate {Tk} by inverting A in two 
steps: first computing Uk = A f 1 (S k ), taking advantage 
of the periodicity of Ai = Ai, and then computing Tk = 


To take advantage of Ai’s periodicity, note that if 8 = Ai(if) 
and u = n + r, with n an integer and r G [0, 1), then 


s = nAi(l) + Ai(r) so that n = and r = A 1 1 (s 

mod Ai(l)). Finally, Af 1 can be efficiently evaluated in 
[0, Ai(l)) by interpolation of a lookup table. 


The input sequence {Sk} can be efficiently generated by 
summing Independent and Identically Distributed (IID) unit- 
parameter exponential random variables Zk, which, in turn, 

may be generated simply as Zk = — log(l — Vk) with Vk ~ 
Uniform[0, 1), (e.g., [25].) 


In summary, this gives the following three step approach to 
event simulation used in the SEXTANT ground testbed and 
end-to-end simulation: 


1. Simulate N standard Poisson process event times 

{Sk}f[ =1 , 

2. Transform these events to phase events Uk = A 1 1 (S k ), 

3. Finally transform to Tk = 0 _1 (T r / c ). 

The photon count increment N should be chosen so that 
the approximation 0(£) ~ 0(£ m ) is valid over the interval 
spanned by the event times {T k } k=1 . If the simulation of a 
continuous observation of a single pulsar needs to be broken 
up to achieve this error condition, the interval breakpoints 
should coincide with a particular event time Tk, to ensure the 
process remains a valid NHPP. 

There are other approaches to simulating NHPPs (see 
e.g., [20]). In particular “accept/reject” methods are simple 
and quite general and may be preferred in some circum- 
stances (time varying particle background rates), but are not 
expected to achieve the efficiency of the method described 
above. 

Filter measurement model 

In the SEXTANT navigation filter, the phase and Doppler 
measurements are blended with models of the spacecraft 
dynamics to update an estimate of the spacecraft state. 

In order to do this, a measurement equation, describing 
how the spacecraft position x and velocity x is related to 
the measurements is needed. For the phase measurement, 
assuming a Geocenter reference observatory (where parallax 
and Solar Shapiro delay are negligible), this is just the phase 
evolution equation 


0(0) = 0o (t - n • x(t)/c). (14) 


7 


previously described in the current section. 

The frequency measurement equation is obtained by differen- 


dating (14), leading to 


4>(t) = 4> 0 {t- A(t)) 


1 - 


n • x(t) 


(15) 


with A (t) := n • x(t)/c. 

The filter also requires the first order partial derivatives which 
are given by 


9y _ 1 f cj)o(t - A(t))ft T 0 \ 

dz c V(1 - A(t))<j> 0 (t - A (t))n T fi 0 {t - A {t))n T J ’ 

(16) 

with y = (0, 0) T and z = (x, x). 


6. GXLT End-to-End Simulation 

The NASA GXLT, developed to support the SEXTANT 
demonstration, is a unique hardware and software test envi- 
ronment that allows for rapid, high-fidelity, end-to-end sim- 
ulation and performance for arbitrary mission concepts and 
evaluation of various spacecraft XNAV scenarios. The GXLT 
leverages several Goddard Space Flight Center (GSFC) Guid- 
ance, Navigation, and Control (GN&C) software tools and 
X-ray detector lab technologies. The overall end-to-end 
simulation architecture of the SEXTANT ground testbed [27] 
is depicted in Figure 6, with the available simulation process 
flows per level indicated by the colored arrows. 

GXLT Simulation Flow 

The GXLT offers three levels of simulation, each with unique 
benefits, where the simulation process flows differ primarily 
in the way measurements are produced. These levels are 
enumerated and described below. 

Level 0: navigation processing only — In this simulation 
mode, the phase and Doppler measurements are directly 
simulated by evaluating the measurement model and adding 
simulated noise. The photon arrival process is not simulated. 
As a standard approach for navigation performance studies, 
this approach is useful for long simulation studies, e.g., deep 
space trajectories, and for tuning navigation filter parameters. 
This flow is indicated by the red arrows in the upper path of 
Figure 6. 

Level 1: software simulated photon arrival process — In this 
mode, the fidelity of measurement generation is increased by 
simulating the photon events in software. The photon arrival 
process is simulated to generate the photon events, which 
are then processed using Eq. (8). This flow is the primary 
development mode for SEXTANT, and is indicated by the 
green arrows in lower path of Figure 6. 

Level 2: real-time, hardware -in- the -loop mode— In this mode, 
which is the subject of Section 7, the MXS is driven in 
real-time to produce X-rays, in the correct energy band, 
with arrival-time statistics that precisely follow the simulation 
truth profile. This dynamically modulated output of the MXS 
stimulates a customized Commercial Off-The-Shelf (COTS) 
detector package, a NICER Engineering Test Unit (ETU) 
flight detector, and eventually NICER flight detector(s). 

Importantly, the modulation of the MXS incorporates the 
orbital dynamics of a simulated spacecraft, so that photon 
events from a static detector emulate a detector following the 


prescribed orbit trajectory Importantly, the MXS is specifi- 
cally modulated so that photon events from a static detector 
appear to originate from a detector following a prescribed 
orbit trajectory. 

The COTS detector package used by the GXLT hardware 
suite includes a commercial Silicon Drift Detector (SDD) 
from Amptek, Inc. and custom timing/readout electronics, 
which log timestamped photon events to a file. These events 
can then be processed by the flight software, and have previ- 
ously been verified to match the statistics of those simulated 
in Level 1 [27]. This flow is indicated by the blue arrows 
along the center-line of Figure 6. 

Simulation procedure 

The simulation process for the flow of each level is described 
in the following. 

Scenario definition — A simulation starts with the specifica- 
tion of a scenario configuration that defines the simulation 
level and length, X-ray optics/detector models, pulsar target 
list and models, target observation times, a truth ephemeris 
file, parameters for the orbit propagator, scheduling options 
and visibility constraints, event simulation options, photon 
processing algorithms, and finally, navigation filter options. 

Truth Trajectory generation — Next, a truth spacecraft 
ephemeris is read or generated based on the scenario options. 
The SEXTANT end-to-end simulation can use an external 
ephemeris file generated by an arbitrary source, or alternately 
run the NASA open source General Mission Analysis Tool 
(GMAT) software to propagate a trajectory on the fly. The 
force models includes a 40x40 EGM96 earth gravity field, 
third-body perturbations for the Sun and Moon, and the 
MSISE90 atmospheric drag model. ISS trajectory data is 
made available to the public by the NASA Johnson Space 
Center (JSC)’s Flight Design and Dynamics Division [28]. 
The available ISS real-time data include: mass (weight), area, 
drag, monthly solar flux and geomagnetic index, maneuver 
information for the current ephemeris, state vectors, two-line 
mean element sets, etc. 

Visibility and Scheduling — Next, the visibility of target pul- 
sars is calculated at each simulation timestep. This process 
considers occultations by the Sun, Earth, and Moon; instru- 
ment constraints that limit achievable azimuth and elevation 
angles; and particle background limits within the South At- 
lantic Anomaly (SAA). 

Then, based on visibility, a pulsar observation schedule is 
generated that attempts to reduce the state residuals by bal- 
ancing pulsar target priorities within pulsar viewing periods. 
Significant observation times per pulsar navigation measure- 
ment, due to signal noise require that observations of each 
pulsar, are grouped to avoid stagnant navigation data. At the 
same time, each pulsar navigation measurement must also 
be timed to arrive when needed to reduce state residuals. 
Instrument hardware slew rate limits are also included in the 
schedule. 

Truth observable generation — At each time step, the detector 
pulse phase model is evaluated using the model at the refer- 
ence observatory and the spacecraft ephemeris. For a geo- 
centric reference observatory, we simply evaluate the phase 
model, as referenced to Universal Coordinated Time (UTC), 
at the reference observatory and account for the Roemer, or 
geometric, delay. For an SSB reference observatory, we trans- 
form the spacecraft state to barycentric coordinates and time 
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Figure 6: SEXTANT ground testbed architecture with three levels of simulation fidelity showing the flow from orbit 
simulation through truth measurement simulation, photon processing and navigation filtering algorithms. 


and evaluate the phase model, as referenced to barycentric 
time, while accounting for parallax, Roemer delay, and Solar 
Shapiro delay. 

Measurement generation — The simulation splits into three 
paths for measurement generation. At Level 0, measurements 
are generated by adding noise to the truth measurements 
consistent with the CRLB, plus a noise component used to 
capture additional noise sources. Individual photon events 
are not generated. 

In Levels 1 and 2, photon events are created that can be input 
into the flight software directly. In Level 1, photon events 
are simulated using the algorithm described in Section 5. 
In Level 2, the photon events are generated in real-time by 
the MXS and timestamped by either a COTS, NICER ETU, 
or NICER flight detector package. Next, in both Levels 1 
and 2, the events are delivered as a sequence of telemetry 
packets to the SEXTANT flight software. Once received, the 
software first pre-processes the events to filter out background 
events that fail a threshold test. Events that pass the pre- 
processor are buffered on a per-pulsar basis. When sufficient 
observation time is accumulated on a particular pulsar, the 
buffer is batch processed as described in Section 5. The 
output is a pulse phase and Doppler measurement suitable for 
processing by the navigation filter. 

Navigation filter — Finally, the pulse phase and Doppler mea- 
surements are used to correct the spacecraft state estimate 
in the navigation filter. The SEXTANT flight software and 
end-to-end simulation use the XNAV-enhanced GEONS Ex- 
tended Kalman Filter for navigation [16, 29]. The ground 
testbed also supports the NASA open source Orbit Deter- 
mination Toolbox (ODTBX) filter software [30]. In either 
package, high fidelity models of orbital dynamics are used 
in conjunction with the measurement model described in 
Section 5 to propagate and update the spacecraft state. The 
propagated state history is stored and fed back to the mea- 
surement generation block to provide the phase predictor 
needed there. 


Figure 7. A key enabling component for the GXLT HWIL 
capability is GSFC’s unique Modulated X-Ray Source. 

Modulated X-ray source 

Most common active X-ray sources use high voltage to ac- 
celerate electrons, emitted from a hot filament, that impinge 
upon a metal target at high velocity. These electrons produce 
X-rays as they decelerate on target impact. Such conven- 
tional electron impact sources are modulated by varying the 
filament temperature or the accelerating potential, and con- 
sequently, suffer a modulation timescale on the order of 1 s. 
Emulating the MSP X-ray environment requires considerably 
faster switching times. GSFC has developed a new kind of 
electron impact source, the Modulated X-Ray Source (MXS) 
seen in Figure 9, which uses an optical or Ultraviolet (UV) 
Light Emitting Diode (LED) to liberate photoelectrons from a 
photocathode [31] that are then accelerated toward a target at 
high potential. Thus, modulating the LED output modulates 
the photoelectron flux, which produces the modulated X- 
rays. The MXS X-ray output modulation is limited only by 
the LED. Commercially available LEDs support switching 
timescales on the order of nanoseconds, which enable the 
MXS to meet our environmental emulation need. 

To simulate the photon arrival process as seen by a detector 
on a moving spacecraft, the MXS input LED is driven with a 
current signal proportional to the count rate function A (</>(£)) 
of the pulsar of interest. The task is achieved by the MXS 
driver. 

Modulated X-ray source driver 

The MXS driver is composed of a control program that 
runs on a PC, the MXS driver digital board, and an LED 
driver circuit, as shown in block diagram form in Figure 1 1 . 
The control program takes as input the pulsar definitions 
and pulsar phase dynamics file, and produces a real-time 
implementation of A (</>(£)) as a precision current signal of 
appropriate level to drive the input of the MXS. This signal is 
precisely aligned to its internal master simulation timer which 
produces the simulation of UTC, t. 

The control program, written in Python [32], runs on a stan- 
dard PC and interfaces to the MXS driver digital board using 
an Ethernet connection. At each simulation time tick, the 
control program feeds target phase commands to the digital 


7. GXLT Hardware-in-the-Loop Mode 

In this section, we describe, in more detail, the Level 2 HWIL 
simulation mode and its hardware components, as shown in 
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Figure 7: SEXTANT Level 2 simulation Hardware-in-the-Loop (HWIL) test configuration. The MXS is driven in real-time to 
produce X-rays with photon events following the simulation truth profile. The X-rays travel through a short gap and impact on 

the X-ray detector and are time-tagged using the FPM/MBU. 



Figure 8: A photograph of the SEXTANT ground testbed. 


board. The digital board is based on a Xilinx ML605 Field- 
Programmable Gate Array (FPGA) card and consists of an 
embedded microprocessor that manages the communication 
with the hardware and custom logic blocks, all implemented 
in the FPGA. The custom logic implements a simulation 
timer, a Numerically Controlled Oscillator (NCO) for gen- 
erating the pulse frequency, a look-up table that stores the 
pulse shape and relative count rates, and an interface to the 
Digital-to- Analog Converter (DAC) that produces the current 
signal which drives a custom LED driver analog circuit that 
ultimately modulates the LED. 

The LED driver circuit, shown in Figure 10, converts the 
differential voltage output of the DAC to a precision current. 
It is composed of a four- stage amplifier circuit that provides 
necessary signal gain, independent controls to adjust the 
signal and background count rates, voltage-to-current conver- 
sion, and safety features to protect the MXS input diode. 

Testing with the pulsar simulator 

The output of the MXS is directed at an X-ray detector. 
A calibration procedure is then executed to ensure that the 
correct absolute count rates are obtained at the detector output 
by adjusting the channel from the MXS to the detector. Typ- 
ically, this will be adjusted so the expected count rate from 
the entire XTI is seen in a single detector. Once absolutely 
calibrated, relative count rate changes upon switching from 
one source to another is faithfully reproduced by the MXS. 



Figure 9: A photograph of the MXS. 


For SEXTANT primary development, the MXS stimulates 
a NICER ETU detectors attached to the XTI timing chain, 
which includes a ETU flight computer that hosts the SEX- 
TANT flight software. This configuration provides for a very 
high-fidelity on-orbit environment emulation. 

In early development, prior to the arrival of the ETU de- 
tector package, a COTS SDD and a custom digital time- 
stamping component, called the Time-Tagging X-Ray Detec- 
tor (XDET), were used to time tag photon detections. The 
XDET consists of a control program running on a PC and 
a digital board board with an Analog-to-Digital Converter 
(ADC) attached. The output of the XDET is a list of photon 
time-of-arrivals stamped with simulated UTC time, which 
are stored in a file and then processed in the end-to-end 
simulation components as shown in Figure 6. 

Figure 13(a) and 13(b) show Level 2 GXLT simulations of the 
Crab pulsar with high and low Signal-to-Noise Ratio (SNR), 
respectively. In these figures, three cycles of the MXS driving 
signal, emulating the rate function A (0(f)) of the Crab pulsar, 
is shown on top of the resulting output of the X-ray detector 
for a low SNR (top (a)) and high SNR (bottom (b)) case. Each 
detected photon results in a impulse shown on the scope trace. 
The rate of events is proportional to the driving signal. 

The MXS driver current signal is shown in blue along with 
the output of the detector, where each individual detected 
photons appear as green vertical lines. Detections occur at 
a rate proportional to the driving signal. Figure 14 shows the 
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Figure 10: Block diagram of the LED driver circuit which is 
used to convert the differential voltage output of the DAC to 
a precision current that drives the input of the MXS. 


MXS driver 



Figure 11: Block diagram of the MXS driver which consists 
of the Control Program, MXS driver digital board, and the 
LED driver circuit. 




(a) 


0 



(b) 

Figure 13: Simulated X-ray events from the GXLT Level 2 
simulation captured on oscilloscope, (a) low SNR case (b) 
high SNR case. 


Folded Light Curve for PSR-J0534+2200 



Figure 14: Simulated X-ray events from the GXLT Level 2 
simulation captured to file and folded using the 
prediction. 


Figure 12: A photograph of a NICER X-ray detector. 


Crab lightcurve clearly emerge as the result of folding the 
resulting events. 


8. Preliminary Simulation Results 

In this section, we present simulation results using the SEX- 
TANT’S ground testbed end-to-end simulation capability, de- 
scribed in Section 6. The simulation described here models 
the SEXTANT primary experiment. 
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background model and integration of background estimation; 
Modeling of additional noise sources including reference 
clock errors, pulsar almanac model mismatch, intrinsic noise 
in the pulsars such as timing noise and glitches, uncertainty 
in the radio to X-ray phase offset, and astrometric errors; 
Enhancements to efficiently handle processing of the Crab 
Pulsar, which is an outlier in the Pulsar catalog, in that its 
flux is 5-6 orders of magnitude higher, and it requires timing 
ephemeris updates much more often than all other SEXTANT 
pulsars. 

The NICER and SEXTANT completed a successful Critical 
Design Review (CDR) in September 2014 and continues to 
make excellent progress. Currently, the first version of the 
NICER flight software has been released for testing on a 
commercial development version of the flight processor, and 
includes the first release of the SEXTANT flight software. 
The second and final release of the SEXTANT XFSW is 
planned for early March of 2015, and will begin testing 
shortly there after. 


Figure 15: Results of the GXLT Level 1 simulation for a 9.5 
day ISS simulation. This plot shows the RSS errors in 
position (top) and velocity (bottom) and their corresponding 
predicted 1 -sigma variances. 


Additionally, all photon events will be telemetered to the 
ground, archived, and made available to the public via NASA 
GSFC’s HE AS ARC. SEXTANT will also use this data to 
support experiments using the GXLT. 


Figure 15, shows results from the Level 1 end-to-end simu- 
lation for a 9.5 day simulation with cyclical measurements 
from four pulsars; B1937+21, J0437-4715, B1821-24 and 
J02 18+4232. These results do not include models of the 
additional error sources (beyond Poisson noise) described in 
Section 5. This simulation is based on fixed observation 
periods of 1800s each broken up over multiple contiguous 
time intervals if needed, taking into account Earth, Moon, 
and Sun occultation constraints. This simulation uses a 
geocenter reference observatory for the phase models (j) q. 
The truth trajectory was generated in GMAT as described in 
Section 6. The GEONS initial state is seeded with GPS-like 
point solution with intentional biases (1km position and lm/s 
velocity each axis) and noise (1km position and lm/s velocity 
lcr each axis) added. The filter uses a JGM-2 30 x 30 gravity 
model, and uses a fit to Harris-Priester drag model, but does 
not estimate the drag coefficient. 

In this simulation (Figure 15), the navigation solution con- 
verges to a level well below the 10km worst-direction objec- 
tive. 


9. Summary and Future Work 

In this paper, we described, in detail, the algorithms and 
models at the core of the SEXTANT technology demonstra- 
tion and the NASA GSFC GXLT that supports SEXTANT 
development, testing, and ground experiments. SEXTANT 
will demonstrate, for the first time, real-time on orbit XNAV 
with a stated primary objective of achieving absolute orbit 
determination to better than 10 km with 2 weeks of measure- 
ments in the highly dynamic LEO of the ISS. While the 
preliminary work contains some optimistic assumptions, the 
early simulation results are very encouraging and indicate that 
the primary objective can be met with margin. 

A number of enhancements to the current baseline XNAV 
system implementation have been clearly identified and are 
in work. These enhancements have been discussed in the 
paper, and are summarized as: Implementation of a variable 
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